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Abstract 

In this work we study the dynamics of the maximum extractable 
entanglement for a system composed of two qubits interacting either 
with two independent thermal baths, a common thermal bath or a 
common squeezed bath. The states with maximum entanglement are 
found applying filtering operations which transform each state to a 
state in Bell diagonal form. We observe a revival of the maximum 
extractable entanglement for common baths. It is also shown that for 
some particular states in two independent baths at zero temperature, 
one can partially recover the initial entanglement at any time. 

1 Introduction 



Entanglement in systems of few particles is one of the most characteristic 
aspects of quantum dynamics. [H [2] . With its counterpart decoherence is a 
key element in quantum computing, quantum cryptography and quantum 
teleportationp]. Quantification of the degree of entanglement pi H] corre- 
sponding to a given quantum state and understanding of how it changes 
due to interactions with other quantum systems or with the environment 
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are of great importance, both from a fundamental point of view or to envi- 
sage possible applications. For systems with many degrees of freedom recent 
developments have stressed approximate methods for the separation of the 
entangled component of a given configuration [21 HI [5], but a complete char- 
acterization has not been achieved. On the contrary for two qubits the works 
of Peres, Horodecki, Hill, and Wootters 0, El El E] established the base for 
a complete discussion of entanglement in terms of algebraic properties of 
the density matrix. In particular, the Peres- Horodecki criterion [6l [7] and 
concurrence as defined by Wootters [9| allow to quantify entanglement for 
arbitrary states. 

When interacting with the environment, entanglement between subsys- 
tems tends to fade away. Interaction may be represented using a bath with 
chosen specific properties. In some cases a complete suppression of entan- 
glement may be observed at finite time, an effect which is referred as entan- 
glement sudden death. [10]. Since in general, interaction with some kind 
of bath is inevitable and one is interested in having maximum entanglement 
at disposal it is important to device strategies to preserve or recuperate en- 
tanglement during the evolution of an open system. Among the strategies 
to minimize the influence of the environment are the use of quantum Zeno 
effect [11], the identification and use of free decoherence spaces [12], and the 
application of quantum error correcting codes in quantum computation [13]. 

In this work we are interested in the use of local operations which allow 
to recuperate at least partially the entanglement which has been lose in the 
interaction with the bath. As is known, there are situations where it is pos- 
sible to improve the degree of entanglement of a pair of systems by means of 
filtering operations consisting in local non unitary operations and classical 
communication (LS0CC)[T3]. Moreover, as we discuss below for each state 
there exists an optimal filtering operation for which the image state is the 
one with maximum entanglement among the accessible sates[T3]. The con- 
currence of this state is what is call the maximum extractable entanglement 
of the original state. In this paper we propose to study the evolution of the 
maximum extractable entanglement for two qubits in contact with a bath 
in order to identify conditions for which local observers using (LSOCC) at 
adequate times may recuperate maximum final entanglement. In the follow- 
ing section we review some of the fundamental concepts discussed in this 
paper. Then in section ([3|) we discuss the explicit form of the optimal fil- 
tering operation. In section (jlj) the dynamics of the maximum extractable 
entanglement for a system composed of two qubits in the presence of two 
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independent thermal baths, a common thermal bath or a common squeezed 
bath is studied. In section the specific case of interaction with the vacuum 
of two independent baths is discussed in some detail. 



2 Concurrence, filtering operations and en- 
tropy 

For a composed system with subsystems A and i? in a pure state, a good 
measure of entanglement is the entropy of either of the two subsystems. The 
entanglement of formation of a mixed state p defined as the average entan- 
glement of the pure states appearing in its decomposition minimized over all 
possible decompositions [15] is also a good measure of the entanglement of 
such system. It is a monotonically increasing function of concurrence intro- 
duced in lU [9] which consequently may be taken as an entanglement measure 
too. For two qubits, with density matrix p , concurrence [9] is calculated in 
terms of the eigenvalues i?2, -R3, -R4 of the related matrix R defined by, 

R = p ay ® ay p* ay (g) ay. (1) 

It is given by 

C = max{{),2^R^-{^^ + + \[Rz + \[Ra)} , (2) 

where 

R^ = max{Ri, R2, R3, R^}- (3) 

Entanglement measured, for example, by concurrence diminishes in general 
by effect of decoherence but may be preserved in some particular situations 
when decoherence free subspaces are allowed. It may also be partially, but 
in general not totally, recovered by applying on the system some specific 
operations. Among them, filtering operations of the form 

{A®B)p{A(g)By 
Tr[{A(g B)p{A® By 

are important because they represent the only way to increase the entangle- 
ment of the bipartite system using local operations and classical communi- 
cation. Under the action of these operations the concurrence transforms as 

^ \det{A)\\det{B)\ 
Tr[A^A®B^Bpy 
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which in particular means that separable and entangled states are mapped 
onto their own kind. 

As we mentioned earlier it is of practical interest to find the state with 
maximum entanglement which can be obtained from the initial state via 
filtering operations. For two qubits, the density matrix can be represented 
in terms of Pauli matrices (a^i t-^ ctq = 12x2, cr*) as 

Ui: <^''^>^^]- (6) 
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It was shown [T3] that the optimal filtering operation maps the initial state on 
a Bell diagonal state. These are states which can be written in the standard 
form 

P=\{l + j:aaf®af'j , (7) 

and define a vector sub-space of three dimensions with coordinates Ci. In this 
subspace, physical states form a tetrahedron whose vertices (Ci,C2,C3) = 
{(— 1, 1, 1), (1, —1, 1), (1, 1, —1), (—1, 1—, 1—)} represent pure Bell states. Sep- 
arable states of this subspace form an octahedron [17]. For matrices written 
in the standard form is easy to show that 

R = p' . (8) 

The eigenvalues of R are {pf, P2, ph pI}, and concurrence in terms of the 
eigenvalues {pi, P2, Ps, P4} of p is, 

C = max{0,2pm- ipi+ p2 + p3 + Pa)} (9) 

where 

Pm = maxipi, p2, P3, Pa} ■ (10) 

In the parametrization ([7]) of the Bell diagonal states the eigenvalues of p 
can be written as 

Pi = ^ (1 + Ci - C2 + C3) (11) 
P2 = ^ (1 - Ci + C2 + C3) (12) 
P3 = ^(1 + ^1 + ^2-^3) (13) 



4 



P, = \{l-Cr-C2-C,) . (14) 

To quantify the degree of mixing of the states during their evolution we 
use von Neuman entropy defined in terms of the eigenvalues pi of the density 
matrix as 

5" = ^ pi In Pi (15) 

i 

For pure states 5 = 0. 



3 Optimal filtering operation 

The practical problem we have to address is given an initial state, to find 
explicitly at any time the maximum entangled state and the optimal filtering 
operation leading to it. The concurrence of this optimal state is call the 
Maximum Extractable Entanglement of the state at this time. There are 
different approaches to compute it which depend on the rank of the initial 
density matrix [TSl [HI [20] ■ Following Leinaas et al [T7j we discuss an exphcit 
procedure to perform this mapping, for a five parameter family of states and 
show how the maximum extractable entanglement evolves in the presence of 
either two independent thermal baths, a common thermal bath or a common 
squeezed bath. Verstraete et al [18] showed that filtering operations (jlD on 
two qubits correspond to Lorentz transformations 



T^J. TU pa 



(16) 



on the real parametrization ([H]) of the density matrix. The Lorentz transfor- 
mations La and Lb are related with A and 5 in (|4]) by 

La = T{A^A*)Ty\det{A)\ 

with T the fixed matrix 



LB=T{B^B*)Ty\det{B)\ (17) 
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(18) 



To identify the optimal filtering operations one should first find [T7] the four- 
vectors m and fi that minimize the function 



F(m, n) 



(19) 
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assuming the normalization condition m^m^j, = n^Uy = 1 and uiq > 0, no > 0. 
Then the optimal Lorentz transformation is chosen so that 



(20) 



In general this transformation does not map directly the initial state to a state 
written in the standard form ([7]). An additional local unitary transformation 
which does not modify the entanglement of the system is necessary to this 
end. Nevertheless for the class of states defined by 0^x4 matrices of the form 
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(21) 



which we consider in this work, the optimal filtering operation leads directly 
to a sate in the standard form. For this class of states using the Lagrange's 
multipliers method the extremal problem (fT9|) may be solved to find the 
Lorentz transformation which diagonalize the real parametrization of the 
density matrix. This is given by, 



(22) 
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where a satisfies both 



and 
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a = ^f^(i±£!L (23) 
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(24) 



with 



p = Vl + a^ . (25) 

Since the Lorentz transformations are not unitary the transformed density 
matrix must be normalized after each step. One then obtains the following 
non zero entries for the final Bell diagonal state: 



13^ + 2af3d + a^c 



(26) 
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Co 



(3^ + 2a(3d + a^c 



(27) 



(28) 



(3^ + 2a(3d + a^c 

Using these expressions and Eqs.(lUl fT^ one obtains the concurrence of the 
optimum state. 



— Free Concurrence 

— Extracted Entanglement 

— Entropy 




Figure 1: Concurrence and maximum extractable entanglement evolution for 
the initial bell state (1,1,-1) which corresponds to the initial state |^I/(0)) = 
H — ) + I — h)) in two independent thermal baths, (n = 0.001) 
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4 Maximum extractable entanglement evolu- 
tion 

In this section we discuss the evolution of the maximum extractable entan- 
glement of a two qubits system in three different situations. We consider 
the two particles interacting with two independent thermal baths, a common 
thermal bath or a squeezed bath. We note that although the final state in 
presence of thermal baths is independent of the initial states, the details of 
the evolution may differ depending of the starting point. Moreover as we 
discuss in detail below for a squeezed bath the steady state also depends on 
the initial condition. 

Since we are interested in configurations with a high entanglement degree 
we choose in each case a pure Bell state as the initial configuration. In what 
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Figure 2: Concurrence and maximum extractable entanglement evolution for 
the initial bell states (1,-1,1) which corresponds to the initial state |^I/(0)) = 
"72'^l~'~~'~)~'~l )) independent thermal baths of n = 0.001 



follows we show the results for the initial states |\E'(0)) = (I H — ) + I — I")) 

and 1^1/(0)) = (I + +) + I ))■ For other possibilities either the initial 

Bell state belongs to the decoherence free subspace and then the evolution is 
trivial or the resulting evolution is similar to discussed cases. 

The master equation for a pair of two-level particles in the presence of 
two independent thermal baths is 

P = 2\^^^ l)(2^aPf^I - (^i(^ap - p(yi(ra) + n{2a\paa - Oaolp - p(7aa\) 
+ {n + l){2abpo-l-o-lahp- po-lab)+n{2alpah-(rhalp- pffbO-l)) (29) 

where (Ta = ® 1, o"f, = 1 ® 0-5 and 7 is the vacuum decay constant. We 
assume that both baths have the same temperature i.e. they have the same 
average number n of thermal photons. 

It is easy to show that evolution with this master equation preserves the 
form of states defined by (12T|) . When Bell states are taken as initial condi- 
tions one obtains for the non zero elements (a, b, c, d) of (1^ a set of equations 
whose solutions can be found either analytically or numerically. In the par- 
ticular case n = the system analytical solutions are discussed in section 
([5]). Meanwhile, using the results of the previous section, we obtain numer- 
ically the evolution of the entanglement and of the maximum extractable 
entanglement. In figures ([1]) and ([2]) we plot the evolution of concurrence 
and maximum extractable entanglement for two different initial states. In 
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the first case we take |\I/(0)) = ^ (| H — ) + | — h)) as initial state and in the 

second the initial state is |\E'(0)) = (| + +) + | ))• 

In this case as in the case with a common bath discussed below sud- 
den death of the entanglement appears sharply. The death time may be 
computed analytically but the result depends on the initial state and is not 
particularly illustrating. In the current case there is a slightly better preser- 
vation of entanglement for the the first initial condition but the entanglement 
evolution is quite similar for both of them. On the contrary the evolution of 
the maximum extractable entanglement is very different. This difference is 
particularly evident for early times. It may be understood as the result of a 
balance between two characteristics of the state: the degree of mixing and 
the entanglement. For a pure state with non vanishing concurrence the ma- 
ximum extractable entanglement is always 1. Correspondingly the filtering 
operations are more efficient on states which are almost pure. On the other 
hand for a mixed state with a very small concurrence the filtering operation 
is not able to improve the degree of entanglement. This is illustrated with 
the use of entropy defined in f|T5|) . In figure ([1]) the entropy decays sufficiently 
fast and the maximum extractable entanglement remains larger until entan- 
glement sudden death when of course it also vanish. For the second case 
(figure ([2])) the entropy for the second initial condition decays too slowly and 
when the state finally is almost pure (entropy almost zero) there is no longer 
sufficient entanglement to be enhanced. 

For a pair of two level particles in the presence of a common thermal bath 
the master equation becomes: 
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P = 2 



{n + l)(2apa^ — aVp — pcrV) + n{2a'^ pa — aa^p — paa'^) (30) 



where a = cTq + (Xf,. The form of the matrix fl2Tl) is also preserved in this 
evolution and as in the previous case one obtains the differential equations 
for to the non zero elements. In figure ([3]) the evolution of concurrence and 
maximum extractable entanglement is displayed for the initial state |\l/(0)) = 

'^l ~' — ^ ~'~ I — figure (jll) for the case in which the initial state 
is 1^(0)) = ^(| + +) + !--)). 

In figure (jl]) we observe that although in this case concurrence always 
decreases there is a sector in which we obtain a revival of the extractable 
entanglement. In other words a decreasing of concurrence does not always 
imply decreasing of maximum extractable entanglement. This is a more 
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Figure 3: Concurrence and maximum extractable entanglement evolution for 
the initial bell state (1,1,-1) which corresponds to the initial state |\l/(0)) = 
(I "I — ) + I — I")) ^ common thermal bath, (n = 0.001) 



dramatic consequence of the competition between entanglement and mixing, 
discussed in the case above. 

Finally we consider evolution of a pair of two level particles in the presence 
of a common squeezed bath. The master equation is 
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P = 2 



n + l)(2crpcr^ — crVp — pcrV) + n{2a^ pa — aa^ p — paa^) 
e'^{2a^pa^ - a^a^p - paV^) - e~'^{2apa - aap - paa) (31) 



where m = yn{n + 1) and ijj are the parameters of the squeezing. 

Once again as in the previous cases, one can verify that the master equa- 
tion preserves the symmetric form fl2T|) of the 6*4x4 matrix. 

One of the most interesting properties of this system is the existence of 
a decoherence free subspace [I2j spanned by the states 



= / , {n\+, +) + -)) , (32) 

102) = ^(|-, +)-!+,-)), (33) 

Each of these states is a stationary state of the dynamics defined by the 
master equation ( 13T1) [21j , [22j , [23j . Moreover for states with no contribution 
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Figure 4: Concurrence and maximum extractable entanglement evolution for 
the initial bell states (1,-1,1) which corresponds to the initial state |^I/(0)) = 
(I + +) + I )) i'^ ^ common thermal bath of n = 0.001 

of the singlet state |02) the system is driven to the pure state |0i). Only for 
initial conditions with non vanishing components along |02) and some other 
direction the final state is a mixed state in the decoherence free subspace |23]. 
This is reflected in the evolution of the system shown in figures ([5]) and 
As one can observe in these figures this system displays two interesting 
effects which are the revival of entanglement and of the maximum extractable 
entanglement. Moreover in the stationary regime almost all the entanglement 
is recovered. This is explained by the fact that for these initial conditions 
there is no component in the direction of |02) and explained above the final 
state is the pure state whose maximum extractable entanglement is 1. 
This is further illustrated by the behavior of the entropy which in each case 
goes to zero confirming that the final state is a pure state. 

5 Maximum extractable entanglement in vac- 
uum n = 

In this section we consider the n = case for two independent baths with 
|\E'(0)) = (I H — ) + I — I-)) as initial state. The non zero components of 
the density matrix satisfies the following equations, 

d = (1 + + 2nd) 
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Figure 5: Concurrence and maximum extractable entanglement evolution for 
the initial bell state (1,1,-1) which corresponds to the initial state |\l/(0)) = 
"75 '^l ~' — ) ~'~ I — ^ common squeezed bath, (n = 0.001) 



a = —7 (a + 2na) 
b = -7 (6 + 2nb) 
c = -27(rf + c + 2nc) 

For n the solutions of these differential equations are, 

-7o(l+2n)t _ 1 

d{t) = 



a{t) 
b{t) 



l + 2n 

,-7o(l+2n)t 
,~7o(l+2n)i 



2e-yo{l+2n)t _ g-27o(l+2n)t _ ]_ 

(1 + 2n)2 







For any finite t when n 

M\t)-{l + c{t)f 
Then from (1231121) and (|2Bll2gl) one obtains 

a ^ 00 Ci ^ 1 C2 ^ 1 C3 



(34) 



(35) 

(36) 
(37) 

(38) 



(39) 



(40) 



The optimal state at any time is equal to the initial Bell state (1, 1, —1). Us- 
ing an infinite boost Eq. fl40l) tells that the complete initial entanglement is 
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Figure 6: Concurrence and maximum extractable entanglement evolution for 
the initial bell states (1,-1,1) which corresponds to the initial state |^I/(0)) = 
"75 '^l ~'~ ~'~ I ^ common squeezed bath of n = 0.001 

recoverable. One can also partially recover the entanglement to any desired 
degree at any finite time using a finite boost. By a numerical analysis one 
verifies that with the boost 7ot ~ y/a one recovers almost all the entangle- 
ment at any finite time t. In figure ([7]) we show the extracted entanglement 
with a = 9. In this case one recovers almost all the entanglement if 7ot ~ 3. 

6 Conclusion 

In this work we studied the evolution of the maximum extractable entangle- 
ment for an open system of two qubits considering three different interactions 
with the environment. For two independent thermal baths at zero tempera- 
ture we show that it is possible to recover almost all the initial entanglement 
using finite operations of local filtering. In the case of a common thermal 
bath we observed during the evolution an increasing of the maximum ex- 
tractable entanglement when entanglement was in fact diminishing. Related 
to this is important to note that the maximum extractable entanglement is a 
property of the state and not of the evolution. In this case what is happen- 
ing is that evolution drove the system to states with less entanglement but 
more extractable entanglement. This suggest as a strategy to manipulate 
efficiently the entanglement to set the conditions of interaction of the system 
with the environment in such a way not to preserve maximum entanglement 
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Figure 7: Concurrence and partial extracted entanglement evolution for the 
initial bell states (1,1,-1) which corresponds to the initial state |^(0)) = 
"75 (I ~' — ) I — ^ common squeezed bath of = 0, a = 9 

but maximum extractable entanglement. 
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